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Recent CDF and D0 measurements of the top quark forward-backward asymmetry at the Tevatron 
present deviations between la and 2a from QCD predictions and a cut M t j > 450 GeV on the 
events increases this deviation to 3.4 a in the CDF measurement. The inclusion of electroweak 
contributions of &{a 2 ) and ff{aa 2 ) enlarges the theoretical prediction from QCD by a factor 
~ 1.2 and diminishes the observed deviations. The calculation method is shortly discussed and 
numerical results are compared to the experimental data. 
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1. Introduction 

In the experimental analyses at the Tevatron, two different definitions of the forward-backward 
asymmetry have been used: 

a (Ay >0)-a(Ay<0) 
o{Ay>0) + o{Ay<0) 

and 

pp _ g(y, >0)-g(y, <0) 

A FB ~ "77 ; n\ , _/ 7777 



A tt _ "l^V ^ V) - uyi±y v. 
A FB — 371 — I ^\ , 777 777T 



where Ay is defined as the difference between the rapidity y t and yj of t and t and the beam axis is 
oriented in the direction of the proton. The values obtained by CDF for the inclusive asymmetry 
[|] are (A§ B = 0.158 ±0.075, A% = 0.150 ±0.055). A% is compatible with the value obtained 
by D0(A% = 0.196 ±0.065) [@]. 

All these values are larger than the Standard Model LO predictions Af B ~ 7%, A P F P B ~ 5% (see 
e.g. [§]) and imposing a cut M ff > 450 GeV, the value obtained by CDF (A% = 0.475 ± 0. 1 14) is 
at 3.4 a from the prediction at this level of accuracy. These results have led to many speculations 
on the presence of new physics and so a thorough discussion of the SM prediction and the corre- 
sponding uncertainty is necessary. At present, the theoretical accuracy is limited by the missing 
calculation of the complete NNLO contribution from QCD to the antisymmetric part of the tt pro- 
duction cross section. Besides the strong interaction, also the electroweak interaction gives rise to 
contributions to the tt forward-backward asymmetry. Although smaller in size, they are not negli- 
gible, and a careful investigation is an essential ingredient for an improved theoretical prediction. 
In the following we briefly summarize our calculation and compare numerical results with experi- 
mental data. This talk is based essentially on [Qj. 

2. Outline of the calculation 

Tree level diagrams of the partonic subprocesses are gluon, photon ad Z s-channel type for 
qq — > tt (Higgs exchange is completely negligible) and s-channel, t-channel and u-channel type 
for gg — > tt. At leading order the production of pairs in pp collisions originates, via the strong 
interaction, from the partonic processes qq — > g — > tt and gg — > tt, which yield the ^(oc s 2 ) of the 
(integrated) cross section, i.e. the denominator of AfB in ( |l.l| ) and (L2). Instead the antisymmetric 



cross section, the numerator of Afb, starts only at &{a^), so the leading term of the asymmetry 
involves one loop corrections to tt pair production. 

Writing the numerator and the denominator of Afg (for either of the definitions ( |1.1[ ) and (1.2)) in 
powers of a s and a we obtain 

N a 2 N + a^Ni + ajaNi + a*N? H — 
Afb = — ~ 



D a 2 D Q + a 2 D ± a s 3 Di ± afaDi + ■ 



=^ +a ^-*w»°> + 4 + g§i + ... (2 „ 

l>o jjq Do a* Uq 

Only some parts of N2 are currently known [Q, ^] and the inclusion of the N\D\/Dq term without 
A?2 would be incomplete, so we have chosen to drop the incomplete &(ccj) part, as done in [j7|]. The 
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inclusion of this term would decrease the asymmetry by about 30%, which indicates the size of the 
NLO QCD term that we dropped 1 . 

The remaining terms include Do coming from the leading &(ccf ) part of the total cross section, 
N\ from the asymmetric part of the NLO QCD corrections to the cross section and No,N\ from 
asymmetric 0(a 2 ), G{a 2 a) parts of the cross section. In the following we show how these terms 
arise and how we (re-)evaluated them (for more details see [0]). 

The squared terms \^ q q^g^ tt t 1 2 and |^# gg ^ w -| 2 yield Do of the LO cross section; the ff(a 2 ) terms 
arise from \^ q q^ r ^ t t + ^qq^z^tt] 2 , which generate a purely-electroweak antisymmetric differ- 
ential cross section, in the parton cms given by 



pe± = 2™ 2 cos e ( i - [rMAA +2K 2 AqAtVqVt 



dcos 8 

1 



Vq = T q 3 -2Q q sm 2 (dw), Aq = T 3 



(2.2) 



4sin 2 (0 H /)cos 2 (%) 

where 6 is the top-quark scattering angle, Q q and Q t are the charges of the parton q and of the top 



and A q , A t and V q , V t are their axial and vectorial couplings to the Z boson. In AfB ( 2. 1 ) this leads to 
the term Nq. The complementary symmetric cross section provides the term Do in the denominator, 
which does not contribute in the order under consideration. Interferences of qq — > y, Z — > ti and 
qq g — » it are zero because of the color structure. Basically for qq — > ti there are also G(a) W- 
mediated £-channel diagrams with q = d,s,b, but they are strongly suppressed by the CKM matrix 
or by parton distributions (q = b). 

The ^(a 3 ) terms that contribute to N arise from four classes of partonic processes: qq — > tt, 
qq — > tig, qg — > tiq and qg — > tiq. In the first case the origin is the interference of QCD one-loop 
boxes and Born amplitudes; the other processes correspond to real-particle emissions. The box 
integrals are free of ultraviolet and collinear divergences, but they involve infrared singularities 
which are cancelled after adding the integrated interference of initial and final state gluon radiation, 
the only asymmetric contribution from qq — > tig at ff(a^). qg — ^ tiq and qg — > tiq yield also 
contributions to AfB, but they are numerically not important [|7p. 

In order to analyze the electroweak G(a 2 a) terms, it is useful to separate the QED contributions 
involving photons from the weak contributions with Z bosons. In the QED sector we obtain the 
&{a 2 a) contributions to Af from these three classes of partonic processes: qq — > ti, qq — >• tig 
and qq — > tiy. The first case is the virtual-photon contribution, which can be obtained from the 
QCD analogue, namely the ^(a 3 ) interference of box and tree-level amplitudes, by substituting 
successively each one of the three internal gluons by a photon, as displayed in Figure [T[ 
In a similar way, also the real-radiation processes qq — > tig and qq —> tiy can be evaluated starting 
from the result obtained for qq — > tig in the QCD case and substituting successively each gluon by 
a photon. 

The antisymmetric &{a 2 a) term from qq — > tig comes from the interference of qq — )• g — )• tig and 
qq — )• 7 — > tig, while in the case of qq — > tiy it comes from the squared amplitude obtained from 
qq —> g — > tiy diagrams. The essential differences between the calculation of the ^"(o; 3 ) and of 



'in a very recent paper [g] the expansion of Afg in powers of a s has been revisited applying the principle of 
maximum conformality. In this case N% and also D2 seems to be numerically negligible, so the term N\D\/Dq can be 
safely included. 
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Figure 1: Different ways of QED-QCD interference at 0{a^a). 

QED ff(aja) terms are only the coupling constants and the appearance of the SU (3) generators in 
the strong vertices. Summing over color in the final state and averaging in the initial state, we find 



that we can relate the QED contribution of the antisymmetric term Ni in (2T) to the &{(X^) QCD 
term Ny for a given quark species qq — > tt(+X) in the following way, 

, . aN? ED 36 a 

Now we consider the weak contribution to Ny. It can be depicted by the same diagrams as for 
qq — > tt and qq — > ttg in the QED case, but with the photon now substituted by a Z boson, involv- 
ing massive box diagrams. The result cannot be expressed immediately in a simple factorized way. 
We performed the explicit calculation including also the contribution from real gluon radiation with 
numerical integration over the hard gluon part. 

Also Z-boson radiation, qq — > ttZ, can contribute at the same order, but it yields only a tiny effect 
of 10~ 5 in AfB and thus may be safely neglected. The same applies to ud — > ttW + as well as to 
Higgs-boson radiation. 

It is important to note that all these partonic subprocesses P1P2 — > tt(+X) can be generated with 
P\ (pi) coming from the first(second) hadron h\(h2) or from /i2 (^i)- Given a kinematic configura- 
tion of pip2 — > tt(+X), if it contributes to o(J t > 0) in the hi Q12) configuration it contributes with 
the same partonic weight also to o(Y t < 0) in the ^(/zi) configuration. So the total contribution to 
Ap P B is non vanishing only if the weights coming from the parton distributions are different, that is 
if: 

fpiM { x \)fp2.h 2 ( x 2) + fpi.h 2 ( x \)fp2M ( x 2) ( 2 -4) 

where f Pi; h (xi) is the parton distribution of the parton pi in the hadron hj. The same argument 
applies also to A^ B . 



At the LHC h\ = ]\2 so AfB, using definitions (1.1) and (1.2), is equal to zero, at Tevatron (2.4) is not 



generally true but it can be used to distinguish which subprocesses can give rise to contribution to 
AfB- Only initial states with at least one of the two p\ and P2 equal to (anti)quark up or (anti)quark 
down can produce an asymmetric contribution. This last statement is completely independent on 
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the assumptions made for the partonic calculation, it relies only on the way proton structure is 
described by partonic distribution functions. 

3. Numerical results 



According to the argument discussed after (2.1), we choose MRST2004QED parton distribu- 
tions for NLO calculations and MRST2001LO for LO, using thereby a s (jl) of MRST2004QED 
also for the evaluation of the cross sections at LO (a similar strategy was employed in [Q]). We 
used the same value pL for the factorization scale and we present the numerical results with three 
different choices for the scale: jJ. = m t /2,m t ,2m t . Other input parameters are taken from [0. 
The various contributions to the asymmetry of either of the two variants A FB and A P F P B are listed 
in Table 0. The ratio of the total ff(aja) + i?(a 2 ) and ^(a $ 3 ) contributions to the numerator N 



(a) A'l, 



(b) A% 
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fl = m t /2 
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fl — 2m t 
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FB 
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fi = mt 


fl = 2mt 


ff{af) 


uii 


7.01% 


6.29% 
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4.66% 


4.19% 


3.78% 


ff{af) 


cid 


1.16% 


1.03% 


0.92% 
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dd 
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QED 
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UU 
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dd 


-0.04% 
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weak 


dd 


-0.03% 


-0.03% 


-0.03% 


G(a 2 ) 


uu 
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0.28% 




(a 2 ) 


uu 


0.11% 


0.14% 


0.17% 


ff(a 2 ) 


dd 


0.02% 


0.03% 


0.03% 


6 


(« 2 ) 


dd 


0.01% 


0.02% 


0.02% 


tot 


PP 


9.72% 


8.93% 


8.31% 




tot 


PP 


6.42% 


5.92% 


5.43% 



Table 1: Different contributions to A" FB andA^g. 



of the asymmetry ([2JJ) gives an illustration of the impact of the electroweak relative to the QCD 
asymmetry. The values obtained for \i = (m t /2,m t ,2m t ) for the two definitions of AfB are 



K EW 



pPP 



V(«l) 

N pp 



N 



PP 



(0.190,0.220,0.254), 
(0.186,0.218,0.243). 



(3.1) 



This shows that the electroweak contribution provides a non-negligible additional part to the QCD- 
based antisymmetric cross section with the same overall sign. Thus it enlarges the Standard Model 
prediction for the asymmetry (the electroweak ff(aja) contribution of uu — > tt to the asymmetry 
is even bigger than the &{ot^) contribution of dd — >■ tt). 



The recent reevaluation of the mixed EW-QCD contribution to Apg in [jlOjj presented values in 
agreement with our results. 

The final result for the two definitions of AfB can be summarized as follows, 



A§ B = (9.7,8.9,8.3)%, 



\pp 

FB 



(6.4,5.9,5.4)%. 



(3.2) 
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Figure || displays the theoretical prediction versus the experimental data. The SM prediction is 
almost inside the experimental la range for Af B and inside the 2a range for A P F P B . It is important 



: 



ml/2 < fi < 2m 
(theory) 



_ 




(a) A% (b) A% 

Figure 2: Theoretical prediction (blue) and CDF data (black=central value, orange=lc, yellow=2a). 

to note that the band indicates the scale variation of the prediction, it does not account for all the 
theoretical uncertainties. For example, as already noted, the &{a^) term in N is missing, and we 
did not include the ^(a s 3 ) part in D. The decrement by about 30% obtained by the inclusion of 
this term can be considered, in a conservative spirit, as an uncertainty from the incomplete NLO 
calculation for the asymmetry. 

(a) A%(M tT > 450 GeV) 



A" 
FB 


ji = mt/2 fx = mt fl = 2m t 


ff{aj) uu 


10.13% 9.10% 8.27% 


ff{aj) dd 


1.44% 1.27% 1.14% 


ff(a 2 a) QED uu 


1.94% 1.95% 1.96% 


ff{a 2 a) QED dd 


-0.14% -0.14% -0.14% 


&(a^a) wea i ( uii 


0.28% 0.28% 0.28% 


0{afa) weak dd 


-0.05% -0.05% -0.05% 


ff(a 2 ) uu 


0.26% 0.33% 0.41% 


ff(a 2 ) dd 


0.03% 0.03% 0.04% 


tot pp 


13.90% 12.77% 11.91% 



Table 2: Different conttibutions to A%(M tt - > 450 GeV). 

We have performed our analysis also applying a cut M t j > 450 GeV to the tt invariant mass. The 
various contributions to the asymmetry A'p B , as discussed above in the case without cuts, are listed 
for M t t > 450 GeV in the Table @. 
The asymmetry with cuts yields 



A" PB (Mti> 450 GeV) = (13.9,12.8,11.9). 



(3.3) 
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(a) A%(M tJ > 450 GeV) 

Figure 3: Theoretical prediction (blue) and CDF data (black=central value, orange=lo\ yellow=2(j). 

A comparison between Table ^a) and Table [j](a) shows that the ratio of the QCD contribution to 
the uu — > tt and old —> tt subprocesses is larger with a M,j > 450 GeV cut, which leads to a slight 
increase of /?f w : 

R tT EW {M t t > 450 GeV) = (0.200,0.232,0.266). (3.4) 

These values of R% w , however, are not enough to improve the situation, indeed the Standard Model 
prediction is at the 3a boundary in case of invariant-mass cut M t j > 450 GeV (see Figure 3 ). 
In Figure || the comparison between theoretical prediction and experimental data from D0 is 
shown. The deviation is larger than in the CDF case (Figure 2(a)), but it is important to stress 
that no statistically significant enhancements have been found by D0 for the region according to 
the cut M t T > 450 GeV. 



4. Conclusions 

The uncertainty of the theoretical prediction for the top quark forward-backward asymmetry 
at the Tevatron is dominated by the incomplete calculation of the contribution from NNLO QCD 
corrections to the antisymmetric cross-section. The electroweak contribution is not negligible and 
increases the LO prediction by a factor ~ 1.2, with differences due to the specific definition of the 
asymmetry and the choice of the renormalization scale. The main part of these corrections is from 
QED origin and it can be derived from the LO contribution multiplied by a simple factor depending 
on the charge of the incoming partons. 

Electroweak corrections cannot explain the enhancement found by CDF including a cut M t - t > 
450 GeV, but they must not be neglected when the deviation is interpreted as presence of new 
physics. 



7 



Top quark AfB-' the elect roweak contribution 



Davide Pagani 



0.35 












0.3 








±2(T 






0.25 




i 




ml/2 <p< 2ml 
(theory) 
















7. 












0.15 


; 


























0.05 






















(a) Ap B 

Figure 4: Theoretical prediction (blue) and D0 data (black=central value, orange=lc, yellow=2a). 
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